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ABSTRACT 
A  procedure  for  deriving  the  variance  of  the  forecast  error  for  Winters' 
Additive  Seasonal  Forecasting  system  is  given.  Both  point  and  cumulative 
T-step  ahead  forecasts  are  dealt  with.  Closed  form  expressions  are  given  in 
the  cases  when  the  model  is  (i)  trend-free  and  (ii)  non-seasonal.  The  effects 
of  renormal ization  of  the  seasonal  factors  is  also  discussed.  The  fact  that 
the  error  variance  for  this  system  can  be  infinite  is  discussed  and  the  rela- 
tionship of  this  property  with  the  stability  of  the  system  indicated.  Some 
recommendations  are  given  about  what  to  do  in  these  circumstances. 

Keywords:  Winters  Additive  Seasonal  Model;  Forecast  Error  Variance; 
Confidence  Intervals;  System  Stability 


All  forecasts  are  wrony.  As  a  result  it  is  rarely  adequate  to  give  a 
forecast  by  itself.  Some  measure  of  the  reliability  or  accuracy  of  the  fore- 
cast is  also  required.  Such  a  measure  is  usually  provided  by  a  confidence 
interval.  This  is  an  interval  within  which  the  future  value  we  are  forecast- 
ing will  fall  with  a  prescribed  probability.  The  philosophy  and  construction 
of  such  intervals  are  well  documented  in  a  number  of  forecasting  texts  and  we 
will  not  consider  them  further  here.  See,  for  example,  Brown  (1962), 
Montgomery  and  Johnson  (1976)  and  Bowerman  and  O'Connell  (1979). 

We  simply  note  that  an  essential  ingredient  in  all  the  expressions  and 
formulae  involved  in  such  intervals  is  the  variance  of  the  forecast  error. 
Consider  the  T-step  ahead  forecast  made  at  time  t  ,  i.e.  the  forecast  of  X 

made  at  time  t  ,  X.(T)  say.  The  corresponding  error  is  e.(T)  =  X.  T   X  (T). 

2  2 

We  shall  denote  its  variance,  var[e  (T)],  by  aT  .  It  is  clear  that  aT  must 

play  a  central  role  in  any  discussion  of  the  usefulness  of  X  (T)  in  predicting 

Xt+T  * 

2 

The  value  of  aT  can  usually  be  computed  on  the  assumption  that  the  under- 
lying model  is  valid.  By  underlying  model  we  mean  the  model  which  the  forecast- 
ing system  assumes  is  generating  the  data.  The  variance  can  certainly  be  derived 
for  the  ARMA  models  of  Box  and  Jenkins,  and  the  interested  reader  is  directed  to 
their  book  (1976).  It  can  also  be  derived  for  all  the  Exponential  Smoothing 
systems  discussed  by  Brown  (1962)  and  others  considered  more  recently  by 
Sweet  (1981) . 

However,  no  results  are  available  for  Winters'  seasonal  systems.  This  is 
particularly  unfortunate  since  these  systems  are  amongst  the  most  commonly  used 
in  practice.  They  are  relatively  simple  to  implement  and  intuitively  appealing. 
Both  systems  are  described  in  detail  in  most  forecasting  texts,  e.g.  Montgomery 
and  Johnson  (1976),  Thornopoulos  (1980),  and  Bowerman  and  O'Connell  (1979).   In 


this  last  work,  approximate  confidence  intervals  are  yiven  for  Winters'  systems. 
As  we  shall  see,  however,  these  intervals  poorly  reflect  the  behaviour  of  the 
error  variance. 

The  Winters'  Multiplicative  seasonal  form  is  a  non-linear  system  and  it  is 
difficult  to  see  how  any  useful  information  about  the  forecast  error  can  be 
obtained  directly.  Approximation  and  simulation  appear  to  be  the  most  sensible 
tools  here.  In  the  case  of  the  additive  model,  however,  we  can  obtain  some 
results. 

The  purpose  of  this  paper  is  to  give  a  procedure  for  deriviny  the  variance 
of  the  forecast  error  for  the  Winters'  Additive  Seasonal  forecasting  system.  We 
deal  with  not  only  point  forecasts  of  future  values,  X,+T  ,  but  also  cumulative 
forecasts,  i.e.  forecasts  of  Y.  T  =  £  X   .  ,  made  at  time  t  . 

Expressions  are  yiven  for  the  variances  of  the  individual  components  of  the 
forecasts,  i.e.  level,  trend  gradient  and  seasonal  factors,  and  for  the  covari- 
ances  between  them.  This  allows  the  construction  of  the  variance  of  any  par- 
ticular forecast  and  hence  of  the  correspondiny  forecast  error.  It  also  enables 
us  to  construct  confidence  intervals  for  the  components  themselves,  and  linear 
combinations  of  them. 

Two  special  cases  are  considered,  viz.  when  there  is  no  trend,  and  when 
there  is  no  seasonality.  This  latter  is  the  well-known  Holt-Winters  non-seasonal 
forecasting  system.   In  both  cases,  closed  expressions  are   obtained  for  the 
variances.  In  the  completely  yeneral  seasonal  case,  the  expressions  derived 
involve  a  number  of  unknowns  which  are  obtained  by  solving  a  set  of  linear 
equations. 

Some  necessary  discussion  is  also  given  about  the  stability  of  these 
forecasting  systems. 


1.  THE  RESULTS 
The  underlying  model  assumed  by  the  Winters'  additive  seasonal  system  is 
given  by 

X  =  m  +  bt  +  sk  +  a  (1) 

where  t  =  rn  +  k  ,  and  n  is  the  length  of  the  season.  The  values 

{s,  :  k  =  1,2,...  ,n}  are  the  additive  seasonal  factors  and  {a.}  is  a  sequence 

of  independent  identically  distributed  random  variates  of  zero  mean  and 

2 
variance  o 

The  standard  form  of  the  corresponding  forecasting  system  is  given  by 

mt  =  ao{h-\J  +  (1-au)(mt-i  +  bt-i} 

bt  =  "l^'V'Vi^  +  ( 1_cxi^bt-i  ^ 

St   =    a2(Xt-mt)    +    (l-a2)St_p    . 

The  T-step  ahead  forecast  is 

Xt(T)  =  mt  +  Tbt  +  St+k.n  (3) 

where  T  =  rn  +  k,  (k  =  1,2,... ,n;  r  _>  0)  . 

The  cumulative  forecast,  i.e.  the  forecast  of 

T  T  . 

Yt  T  =  I     Xt+i   is  given  by     Yt  T  =  I     Xt(i)  ' 
t,i   i=1  L+i  t,i   1=1 

Our  derivation  of  the  forecast  error  variances  parallels  the  development 
given  for  Exponential  Smoothing  models  by  Brown  (1962).  By  definition,  X  (T) 
depends  on  only  past  values  of  a   and  so  is  independent  of  Xt+T  .  Thus, 
ol   =  a~   +  var[X.(T)J  .  Further,  X  (T)  can  be  expressed  as  a  linear  combinatior 

I  L  ^ 

of  these  past  values  of  a   and  so  its  variance  is  proportional  to  a  .  Thus, 


Oj  =   a   [1  +  V(T)]  .  When  dealiny  with  cumulative  forecasts,  the  same  aryument 
holds.  The  error  now  is  E.  T  =  Y.  T  -  Y.  T  and  we  may  show  it  has  variance 

of  the  form  a?'U   +  V£(T)]  . 

2 
In  practice,  a   is  unknown.   It  may  be  estimated  from  the  data,  however, 

2  2 

using  the  fact  that  a  [1  +  V(l)]  =  a,  ,  the  variance  of  the  one  step  ahead 

forecast  error.  Thus,  we  can  estimate  a,~  directly  from  the  forecast  errors, 

and  knowledge  of  V(l)  yields  an  estimate  of  a     . 

We  derive  V(T)  and  VV(T)   from  the  variance-covariance  matrix  of  the 

estimates  of  the  individual  components  of  the  forecast  i.e.  m  ,  b  and 

S^    ,  k  =  1 ,2 , . . .  ,n  .  These  variances  are   of  some  interest  in  their  own 
t+k-n 

right  for  the  construction  of  confidence  intervals  for  the  components,  and  we 
display  them  below.  Details  of  the  derivations  are   given  in  Appendix  1. 

oZm  =   Var(m  )  =  a(j{a*[nd  +  2   £  (n-i)d.J  +  2(  1-^)  (dy-d  )  }a2 

l-l 

2  2  2       2 

ab  =  Var(bt)  =  2aua1(d()-dn)a 

a$  =  Var(St+R_n)  =  2^(1-0^)  (dQ-d1)a  ,  k  =  1,2,... ,n 

°mb  =  Cov(mt'bt)  =  a0al(2"al)(d0-dn)a2 

%k  =  Cov(bt»St+k-n)  =  a0<xla2(1'a0){dnrk_dk  +  dk-l  "  Vk+1^'  k  =  *»2>-<-»n  • 
amk  =  Cov(mt»St+k-n)  ~~   a0a2(1-a0)C(1-al)(dn-k-dk)  +  (dk-rdn-k+l )]a? '  k  =  l'2»'v.n 
aij  =  aji  -   Cov(St+i-n'St+j-n)  =  a2(1-a0)^(2dj-i  -  dj-i-l  -  dj-i+l)o2>  l  1  j  <  J  1  n 


Before  discussing  the  derivation  of  the  values  {d  ,d..,...,d  }  we  note  the  forms 
of  V(T)  and  VC(T)  . 


Now,   V(T)o     =  var(m.    +  TX   +  S.,        ) 

t  t         t+k-n 


2         2  2         2 

=   °m  +   T   o,     +    a     +   2(a  ,     +   Ta  .    +   Ta  , 
m  b         s  mk  mb  bk' 


Thus, 

V(T)    ■   2aQ^(d0-dn)T2 


+  ZVlC<\>(2-"l)(d<rdn)   +  «2{l-«0){6n-k  ~  dk  +  dk-l   "  dn-k+l)]T 


+  2a0a2(l-a0)[(l-a1)(dn.k-dk)   +   (d^-d  )] 


n-1 


+   2a^(l-a1)(d0-dn)    +   2a2(l-au)2(d0-d1)    +   a^[nd0   +   2      J      (n-i)d-] 

i  =  l 

In  the  same  way,   using  the  fact  that 


Yt,T=  Tmt  +  1  T<T+1>bt  +  X   St+k-n 

k=l 

we  can  show  that  when  T  =  pn  +  q,p>_0,l<q<n: 


V£(T)    =  [T2a2  +  T2(T+l)amb  +  ±  T2(T+l)2a2]/o2 

+   2a2(l-a0)2[(p+l)2(d0-dq)    +   P2(VVq>    +   P<P+1>  <  VWq'V  ] 
+   ^ajaja-agJCUp+lM^-dJ    +    (d        -d    )]T(T+1) 


*   2a0a2(l-a0)[p(2-ai)(d0-dn)    ♦    (d„-d   -d,+d        ) 


+  «1     I     <Vdn-lc»T 
k=l 

Now  the  numbers  dr.,d1,...,d   are  necessary  to  evaluate  the  variances. 

0  1     n  J 


(5) 


These  (n+1)  values  are  derived  by  means  of  an  algorithm  which  is  given  and 
discussed  in  Appendix  2.  It  is  a  fairly  straightforward  algorithm  due  to 
Wilson  (1979).  It  has  the  property  that  it  also  checks  the  stability  of  the 
forecasting  system.  As  we  shall  see  in  Section  3  this  is  most  useful  for  this 


forecasting  system.  Before  discussing  any  of  the  technical  aspects  ot  the 
procedure  described  here,  we  present  some  particular  cases  which  have  closed- 
form  solutions. 


2.  SPECIAL  CASES 
2.1.  Trend-free  additive  seasonal  model. 

In  this  case  the  model  assumed  is  a  constant  level  plus  addi- 
tive seasonals.  The  forecasting  equations  are  obtained  by  removing  the  second 
equation  from  the  original  three  (2),  and  dropping  b  .  from  the  first. 
Here  the  values  d  ,d.,...,d   can  be  derived  algebraically.  The  van'ance- 
covariance  matrix  of  the  components  is  now  obtained  from 

2   -  n  2h 

am  "  na0h  » 

as  =  o^d-oiQKln-l)^  +  a2(l-aQ)]h  , 

amk  =  aoa2(1_Vn  '       l   -  k  -  n  ' 

o .  ■   =   a  ■  ■    =   -    a  ,     ,  l<i<j<n, 

lj         Ji  mk    '  -  J  -       ' 

where     h  =  a  VL^-aQ-ap+aQa^)  (naQ+ao-oiQou)]   . 
Thus, 

1   -    (l-an)(l-a?) 
^l^l-a^l-a')    '  <«> 

It  is  interesting  to  note  that  for  this  model  V(T)  is  independent  of  T  . 
For  the  cumulative  forecast  error,  we  find  that  for  T  =  pn  +  q,  (p  >_  0, 

1  <  q  <  n),  VE(T)  =  PT2  +  ^(l-a^^pT  +  P^  +  ^ 

2 
where  P  =  aQ^aQ  +  apd-ajjh/a^ 

2.2.  Non-seasonal  Holt-Winters'  system. 

The  appropriate  equations  are  obtained  by  removing  the  third  equation 

from  (2)  and  S.    from  the  first.  This  is  a  well-known  non-seasonal  fore- 
t-n 

casting  system  and  often  used  for  linear  trend  models.  It  is  a  direct  compet- 
itor of  Brown's  second  order,  or  double,  exponential  smoothing  which  we  will 
consider  shortly.  The  underlying  model  is  assumed  to  be  a  simple  linear  trend 

i.e.  s   is  removed  from  (1). 
k 


We  obtain 

°m  =  (2al  +  2%  "   3a(jal)hl  • 
ab   =  ^«ua1h1  , 

0mb  ~~   «oal(2~al)hl  ' 
wnere  h   =  a  /  (4-2a.-oua,  )  . 

We  can  deduce  in  the  usual  way  that 

V(T)  =  [2aQa2T2  +  2(^^(2-^)1  +  ( 2^+2^-3^)  ^/a2         (7) 

and 

VE(T)  =  CaQa2T2  +  4^0^  +  ^a^a^a^-a^2)  l^h^a2 


2.3.  Brown's  Double  Exponential  Smoothing. 

It  is  well-known  that  this  system  is  equivalent  to  the  Holt-Winters' 

2 
non-seasonal  form  above  if  we  choose  cu  =  1-3"  and  a-,  =  ( 1- 3) / ( 1+3)  where 

3  =  1-a  and  a     is  the  smoothing  constant.  Thus,  the  corresponding  results 

can  be  obtained  by  making  these  substitutions  in  the  equations  of  Section  2.2. 

2   2 

above.  The  results  for  a   ,  a,  ,  a   .   and  V(T)  may  be  found  in  Brown  (1962) 

m   b   mb  J 

and  V,-(T)   is  given  in  Bowerman  and  O'Connell  (1979). 


2.4.  Continuously  re-normalized  seasonal  factors. 

In  the  general  additive  seasonal  model  given  by  (1)  it  is  usually 
n 
assumed  that  Y  s,  =  0  .  This  is  done  to  ensure  some  measure  of  independence 

k=l  K 
between  the  level  of  the  data  and  the  seasonal  pattern.  It  is  always  recom- 
mended that  when  the  corresponding  forecasting  system,  as  given  by  (2),  is 
implemented,  the  seasonal  factors  sum  to  zero  initially.  What  is  to  be  done 


thereafter  is  not  so  clear.  Because  of  the  revision  equation  for  the  seasonal 
factors  they  will  no  longer  sum  to  zero  after  the  first  observation.  We  can 
renormalize  the  seasonal  factors  at  any  time  by  subtracting  from  each  the 
average  of  the  set,  i.e.  the  most  recent  n  .  Whether  we  should  do  so  or  not 
is  not  clear.  Some  authors  seem  to  recommend  against  it,  e.g.  Bowerman  and 
O'Connell  (1979);  some  regard  it  as  an  optional  modification,  e.g.  Montgomery 
and  Johnson  (1976);  some  recommend  renormal ization  once  per  season,  e.g. 
Chatfield  (1978);  and  some  recommend  continuous  renormal ization,  i.e.  after 
every  revision  of  seasonal  factors,  e.g.  Thomopoulos  (1980). 

Our  only  interest  here  is  in  how  such  a  procedure  affects  the  forecast 
error  variance.  In  the  cases  of  occasional  or  purely  seasonal  renormal ization 
the  situation  is  very  complex  and  we  have  nothing  to  say,  except  that  the 
effect  on  the  error  variance  appears  to  be  small.  In  the  case  of  continuous 
renormalization,  however,  the  following  result  may  be  applied.  An  outline  of 
the  proof  is  given  in  Appendix  1.  Continuous  renormalization  of  the  general 
forecasting  system  given  by  (2)  yields  exactly  the  same  forecasts  (and  so 
forecast  errors)  as  running  the  system  without  renormalization  but  replacing 
a0,a-,,a2  by  an  .at  cu  respectively,  where  o^  =  o^  -  a^l-a^/n  , 
a*  =  olqolJql*    ,  and  a*  =  cu/(l+cu/n)  .  As  a  consequence,  the  appropriate 
covariance  matrix  elements  and  values  of  V(T)  and  V  (T)  can  be  obtained 

JL.      JL      Jb 

as  in  Section  1  above  by  replacing  cu,a,  ,ou  by  cu,a,  ,cu  respectively.  The 
same  holds  true  for  the  trend-free  version  discussed  in  Section  2.1. 


3.  TECHNICAL  CONSIDERATIONS 
3.1.  Infinite  Error  Variance  and  Stability. 

In  smoothing  systems  in  general  the  variance  of  the  forecast  error 
increases  as  the  values  of  the  smoothing  constants  increase.  For  example,  in 
the  case  of  the  simple  exponentially  weighted  moving  average  (SEWMA), 
X.  =  aX.  +  (l-a)X.  -.  ,  the  variance  of  the  forecast  error  is  proportional  to 
2/(2-a)  which  clearly  increases  with  a  in  (0,1).  The  need  for  higher  values 
of  a  in  practice  reflects  the  fact  that  the  underlying  level  is  changing 
rapidly.  Consequently,  a  more  responsive  forecast  is  needed.  Moreover,  the 
inherent  instability  of  the  underlying  model  is  reflected  in  the  increase  in 
forecast  error  variance.  This  variance  is  finite  while  a  remains  in  the 
stability  region  of  the  system,  i.e.  (0,2)  for  the  SEWMA.  Note  that  the  set 
of  values  from  which  a  is  usually  chosen  is  a  subset  of  the  stability  region 
of  the  system.  Thus,  a  stable  system  always  results  for  the  SEWMA.  This  is 
also  true  for  all  the  General  Exponential  Smoothing  models.  However,  it  is 
not  true  for  the  seasonal  system  under  consideration  here. 

It  is  a  somewhat  surprising  and  problematic  fact  that  there  are  choices 
of  the  smoothing  constants  an ,  a,,  a?  lying  in  the  usual  range,  (0,1),  which 
yield  an  infinite  variance  for  the  forecast  error.  The  algorithm  given  in 
Appendix  2  checks  for  this  possibility  which  indicates  that  the  forecasting 
system  is  unstable.   If  such  a  situation  arises  clearly  no  meaningful  confi- 
dence intervals  can  be  constructed.  More  importantly,  however,  we  must  decide 
how  to  interpret  this  knowledge  of  the  system's  instability.  The  concept  of 
stability  for  a  system  of  difference  equations  is  an  important  one  but  rcirely 
discussed  in  the  context  of  forecasting  systems.  Two  useful  exceptions  to 
this  are  the  papers  by  McClain  (1974)  and  Brenner  et  al ,  (1968). 
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In  essence,  a  forecasting  system  such  as  (2)  above  is  stable  if  the 
influence  of  earlier  observations  decreases  with  the  passage  of  time.  Thus, 
the  forecasts,  or  (equivalently)  the  forecast  errors,  are  influenced  more  by 
recent  observations  than  by  those  in  the  past.  It  is  worth  notiny  that  this 
is  also  the  essential  philosophy  of  exponential  smoothiny  systems,  and  indeed, 
most  forecastiny  systems. 

The  converse  of  this  is  that  in  an  unstable  system  past  observations  have 
a  constant  or  even  yrowiny  influence  on  future  forecasts.  As  an  illustration, 
consider  the  SEWMA  ayain.  The  forecast  X.   can  be  written  as  a  weiyhted 
averaye  of  all  past  observations.  At  time  t  ,  the  weiyht  yiven  to  X  .   is 

L  —  K 
If 

a(l-a)   .  Clearly,  if  a  lies  outside  the  stability  reyion  (0,2)  this  weiyht 
increases  with  k  so  that  data  in  the  most  remote  part  have  yreatest  influ- 
ence upon  the  forecast.  Equally,  if  a  =  2,  all  observations,  however  distant 
in  time,  make  the  same  contribution  to  the  forecast. 

The  inescapable  conclusion  is  that  it  would  be  extremely  unwise  to  select 
smoothiny  constants  which  do  not  lie  in  the  stability  reyion  of  the  system. 
On  the  one  hand,  we  can  tolerate  a  hiyh  (but  finite)  forecast  error  variance 
because  this  represents  a  trade-off  between  accuracy  and  robustness.  In  an 
effort  to  predict  a  model  whose  parameters  are  chanyiny  rapidly  in  time  we 
may  require  a  more  responsive  system.  The  cost  of  this  is  a  correspondi nyly 
hiyher  error  variance.  On  the  other  hand,  we  can  not  tolerate  an  infinite 
error  variance  because  this  indicates  an  unstable  forecastiny  system.  Such  a 
system  violates  exactly  those  assumptions  which  we  hold  most  important  to  the 
yeneration  of  our  forecasts. 

3.2.  Practical  consequences. 

It  is  obvious  from  the  forms  of  V(T)   in  the  two  special  cases  (6)  and 
(7)  that  the  usual  ranye  of  values  of  smoothiny  constants,  i.e.  (0,1),  lies 
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within  the  stability  regions  for  these  systems.  The  problem  arises  for  only 
the  general  seasonal  model.  The  actual  stability  region  is  difficult  to  es- 
tablish in  general  here  for  it  depends  not  only  on  ou,  ««  and  ^  but  also 
on  n  ,  the  length  of  the  season.  However,  Gardner  (1984)  reports  that  Sweet 
has  demonstrated  numerically  that  the  usual  range  for  the  three  smoothing 
constants  lies  within  the  stability  region  for  n  up  to  four.  Kor  seasons 
longer  than  four  periods  it  is  no  longer  true:  a  result  of  importance  for 
weekly  and  monthly  data. 

The  values  of  the  smoothing  constants  are  established  in  one  of  two  ways 
in  general.  They  may  be  selected  intuitively  by  appealing  to  the  ideas  of 
required  speed  of  response  and  constancy  of  the  underlying  model.  Alterna- 
tively, they  may  be  chosen  as  giving  the  best  fit  with  respect  to  some  cri- 
terion such  as  least  squares.  However  the  values  are  chosen,  we  should  be 
very  concerned  if  they  do  not  lie  in  the  systems's  stability  region. 

If  we  have  chosen  them  for  rapid  response  we  may  allow  their  use  for  a 
brief  period,  perhaps  when  we  initialize  the  forecasting  system.  As  noted 
above,  however,  if  we  do  not  replace  them  with  values  in  the  stability  region 
then  the  longer  we  forecast  the  more  influential  become  the  very  first 
observations.  If  some  best  fitting  criterion  selects  'unstable'  smoothing 
constants,  i.e.  ones  leading  to  an  infinite  error  variance,  then  a  very  real 
possibility  is  that  the  model  is  wrong,  at  least  for  part  of  the  data.  In 
particular,  it  may  well  be  that  a  multiplicative  seasonal  model  is  called 
for. 

We  may  decide  to  retain  the  model  but  use  smoothing  constants  in  the 
stability  region.  These  will  be  obtained  by  decreasing  one  or  more  of 
ol,  ,  a,  ,  a?  .   In  general,  if  a(),  a,,  a-  do  not  lie  in  the  stability  region 
it  will  be  because  one  or  more  of  them  is  too  large. 
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4.   SUMMARY 
The  purpose  of  this  paper  has  been  to  give  a  procedure  for  deriving  the 
variance  of  the  forecast  error  for  Winters'  Additive  Seasonal  System.  In 
summary  that  procedure  is  as  follows: 

(i)   If  the  model  has  no  trend  or  is  non-seasonal  the  results  are  given 
in  Section  2.1  and  2.2  respectively  in  closed  form, 
(ii)  From  the  chosen  smoothing  constants  <*n,  a,  and  o.?     the 

coefficients  W  ,  W  ,...,  W     are  derived  as  in  Appendix  1. 
(iii)  The  W.  values  are  used  to  start  the  algorithm  of  Appendix  2  and 
yield  d  ,d,,...,d  .  The  algorithm  simultaneously  checks  the 
stability  of  the  forecasting  system  for  the  chosen  values 
cxq,  a-,,  cto  .  If  n  £  4  the  system  will  be  stable, 
(iv)  The  values  of  d.  are  substituted  into  the  appropriate  expressions 
for  V(T)  and  V  (T)  . 

Note  that  at  step  (iv)  we  are  able  to  evaluate  the  variance  of  the  errors  for 
T-step  ahead  forecasts  and  T-step  ahead  cumulative  forecasts.  Confidence 
intervals  for  these  future  values  can  be  derived  from  these  variances. 

A  discussion  of  how  to  calculate  the  corresponding  results  when  the 
seasonal  factors  are  renormalized  after  each  observation  is  given.  We  have 
also  discussed  the  (real)  possibility  that  the  error  variance  may  be  infinite 
and  indicated  its  relationship  with  the  stability  of  the  forecasting  system. 
Some  recommendations  are  made  about  the  interpretation  of  system  instability. 
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APPENDIX  1 
The  results  of  this  paper  are  obtained  using  the  fact  that  the  forecast- 
ing system  given  by  (2)  has  an  equivalent  ARIMA  process.  The  nature  of  this 
equivalence  is  discussed  in  McKenzie  (1984)  and  the  equivalent  ARIMA  for  (2) 

is  given  in  McKenzie  (1976).  It  has  the  form  ( 1-B ) ( 1-Bn ) Xf  =  W(B)et  ,  where 

k  n+1    k 

B  is  the  backshift  operator  defined  by  B  X  =  Xt_k  and  W(B)  =  1  -  J  W  B 

k  =  l 

where     kL    =   1   -   cxQ -   a^    ,   Wn   =    1   -    OgC^   -    o^d-o^)    ,   Wn  +  1   =   -    ( l-aQ)  ( 1-a.,)    , 

and  W  =  -  a.,  a,   for  k  =  2,3,...,n-l  . 

00 

Suppose  now  that  W(B)6(B)  =  1  where  6(B)  =  J  6,Bk  .  Usiny 

k=0  K 
equations  (2),  it  can  be  shown  that,  for  the  purposes  of  the  variance  calcula- 
tion, the  components  of  of  the  forecast  can  be  expressed  as  infinite  moving 
averayes  in  {a  }  .  Thus: 


m 


t   =   a0{l+    ^    U.    -    (l-a1)ci.1]B1}at    , 
bt  =   WX   6iBi   +    i     (6i-5i-n)Bl]at    • 


i=0 


i=n 


St  =  ^(I-Oq)!!  +  J     (6i-6i_1)B1]at  , 

min(k,n-l) 
where  5.  =  £  6,  .  .  Further,  the  corresponding  representation  of  S.  . 

n  —  k 
can  be  obtained  by  writing  it  as  B   S.  ,  (k  =  l,2,.,.,n)  . 


Since  {a  }  are   independent  random  variables  the  variances  and  covariances 

can  be  obtained  directly  from  these  moving  average  representations.  Defining 
00 

di<  =  I     ^i^i  +  k  '  ( ^  =  u>l>*'*>n)>  yields  the  expressions  yiven  in  the  paper. 

1=0 
The  alyebra  is  tedious  and  hardly  illuminatiny  and  so  is  omitted.  It  can  be 

obtained  upon  request  from  the  author. 
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As  regards  the  renormal ization  procedure,  recall  the  equivalent  ARIMA  for 
(2)  yiven  above.  If  we  introduce  continuous  renormal ization  as  discussed  we 
find  there  is  still  an  equivalent  ARIMA  and  it  has  the  same  form.  Now, 
however,  W(B)  has  cu,  a,,  cu  replaced  by  a. ,  a*  at    .  This  ARIMA  is  equiv- 
alent to  the  system  (2)  with  the  starred  smoothing  constants.  Hence,  the 
result.  Initially,  the  result  is  surprising.  However,  note  that  we  are 
really  dealiny  with  two  different  decompositions  of  the  seasonal  factor.  One 
form  is  (m.+S. +. _  )  as  it  appears  in  (2)  and  the  other  is  (mt+^t+k-n^  wnere 

S*     is  normalized  to  sum  to  zero  and  m*  is  the  correspondingly  adjusted 
t+k-n  t 

level . 


lb 


APPENDIX  2 
To  obtain  the  values  {d,  }  note  that,  by  definition,  they  are  the  variance 
(d  )  and  the  first  n  autocovariances  (d- ,d?, , , .  ,d  )  of  the  autoreyressive 
process  of  order  (n+1)  given  by  W(B)Z  =  e   ,  where  {e  }  are  independent 
random  variables  of  zero  mean  and  unit  variance.  Thus,  the  sequence  {d  }  may 
be  obtained  by  solving  a  suitable  set  of  the  Yule-Walker  equations  for  this 
process.  This  procedure  is  discussed  in  McLeod  (1975,  1977). 

From  our  point  of  view,  however,  a  much  superior  approach  is  presented  by 
Wilson  (1979).  The  stability  of  the  forecasting  system  corresponds  to  the 
stationarity  of  the  autoregressive  process  {Z  }  and  can  be  tested  routinely 
within  the  procedure.  The  alyorithm  is  as  follows: 

(i)  define  Wp+1  k  =  Wk  ,  k  =  1,2,..., n+1;  and  tR+1  =  1  . 
(ii)  apply  the  following  equations  in  the  given  order  for 
k  =  n+l,n,.. .  ,2 

Dk  =  !  -  Wk,k 

if     D.    —  0    »  tne  system  is   unstable;    stop  now. 

wk-l,i  -  (wk,i  +  Wk,kWk,k-i)/Dk'  <1  -  1.2..-.k-« 

Vl  =  VDk 
(iii)     Di   =   !  "   wisi 

if  D.  <_  0  ,  the  system  is  unstable;  stop  now. 

d„  =  t,/D,  ,  and  the  system  is  stable,  and  the  error  variance  finite. 

(iv)  dj  =  WU1  dQ 

k-1 


lk  =    J=l   Wk-l,i  dk-i  +  Wk,k  \    '  k  =  2,3.  —  ." 


Note  that  if  n  <_  4  the  system  is  stable  for  the  usual  choice  of  parameters 
i.e.  a-  e  (0,1),  i  =  U,l,2,  and  so  we  need  not  test  D.  . 
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